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Abstract 

The fluctuation-dissipation theorem together with the exact density response 
spectrum for ideal quantum gases has been utilized to yield a new expression for 
the static structure factor, which we use to derive exact analytical expressions for 
the temperature-dependent pair distribution function g(r) of the ideal gases. The 
plots of bosonic and fermionic g(r) display "Bose pile" and "Fermi hole" typically 
akin to bunching and antibunching as observed experimentally for ultracold atomic 
gases. The behavior of spin-scaled pair correlation for fermions is almost featureless 
but bosons show a rich structure including long-range correlations near Tc. The 
coherent state at r=0 shows no correlation at all, just like single-mode lasers. The 
depicted decreasing trend in correlation with decrease in temperature for T < Tc 
should be observable in accurate experiments. 
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Introduction: The surge in the study of various properties of ultracold atomic gases 
has prompted search for atomic analog of the Hanbury-Brown-Twiss (HBT) effect [1] 
by various research groups as reported in Refs. [21 El H] and references therein. The 
suppression of density fluctuations, a signature of the Pauli exclusion principle at work in 
real space and thereby antibunching, has been demonstrated in jSJiSj. On the other hand, 
Jeltes et al. |1] have compared results of the two-particle correlations of a polarized, but 
not Bose-condensed, sample of ultracold ^He* atoms with those of polarized ^He* atoms. 
The experimental conditions in [1] were such that the gases could be treated almost 
ideal. Hence bunching for bosons and antibunching for fermions at small interatomic 
separations have been attributed to purely quantum effects associated with the exchange 
symmetries of wavefunctions of indistinguishable particles. Also, the measurement of 
correlations has been reported |5] both above and below the Bose-Einstein condensation 
(BEC) temperature in atomic ^He*. 

An elegant form for the correlation function v{r) of the density fluctuations in ideal 
quantum gases (IQG) has been derived by Landau and Lifshitz [i6j. Although many qual- 
itative features and limiting expressions for z/(r), which is related to the pair distribution 
function (?(r), have been discussed by them, quantitative descriptions would require avail- 
ability of general analytical forms. References [71 El E] have discussed expressions for 
g^^{r) of an ideal Bose gas (IBG) in one or the other range in the temperature domain 
0<T<T+ whereas Lee and Long [10] have given g^^{r) of an ideal electron gas at T=0. 
Reference [IT] has discussed model analytic expression for the unpolarized homogeneous 



electron gas in solids. However, these expressions cannot be utilized to get comprehensive 
theoretical values to compare with the observed HBT effect reported in |1] and [5]. The 
main purpose of this Brief Report is to fill up the gap by deriving exact analytical expres- 
sions for g{r) of IQGs. The improvements on some asymptotic results, both for fermions 
and bosons, available in the literature are also discussed. 

A unified approach is presented for evaluation of temperature-dependent g{r) for ideal 
Bose-Einstein (BE), Fermi-Dirac (FD) and Maxwell-Boltzmann (MB) gases wherein the 
unification has been achieved using polylogarithm [121 1131 [IH IISI [161 HI]- Our method 
starts from the expression for the density response function derived in [17] using the 
method of second quantization, applies the fluctuation-dissipation theorem to get a gen- 
eral new expression for the static structure factor, introduces a function which is related 
to the one-particle density matrix, and ultimately gets the general form for g{r) valid 
for all quantum gases at all temperatures. The computed values of g^^{r) and g^^{r) 
using our analytical expressions are depicted graphically as a function of r at various 
temperatures showing at small-r appearance of bump and dip which we term as "Bose 
pile" and "Fermi hole". The plots are further compared with the experimental results 
[H [5] for ultracold atomic gases. 

Basic Expressions: The pair-distribution function g{r) of a uniform one-component fluid 
consisting of particles in volume V is defined by the thermal average of an operator 
that counts pairs of particles located distance r apart, divided by the square of number 
density. It is related to the static structure factor S{q) of the fluid by the spatial Fourier 
transform: 

n[9ir)-l] = We'^- [S{q)-1] (1) 
^ q 

with n = N/V denoting the number density. Taking due account of the fact that the 
operator of the total particle-number, = A^q=Oj is a constant of motion, the fluctuation- 
dissipation theorem for a uniform system, x"(g, w) = {mT/h){l — e~^^'^)S{q,uj), can 
be solved for the van Hove function S{q,uj) = (fi/nir) (1 — (5q^o) x"{(li^)/ (l ~ e"''^'^^ + 
(^{SN)'^ /N^j 6fifi6{u)), where (...) represents averaging in the grand canonical ensemble 
(GCE). The expression S{q) = /_f^da; S{q,u) then yields 

S{q) = (1 - <5q,o) ^ fdu coth (^) x"iq, ^) + '^q.o^^^^- (2) 



Upon inserting Eq. (I2|) in Eq. ([T|), the q-sum on the right-hand side separates into two 
parts: the first part contains the summation with the restriction q 7^ while the second 



which vanishes in the thermo- 



'n ( {5NlNf ) - l/V 

dynamic limit for all IQGs except for an IBG at T < T^. The pathological aspect of 
GCE for the condensate fluctuations has been ameliorated by replacing ( [(5(ao'3o)]^ 



results in an additive constant, 



GCE 

by [(5(agao)]^ ^ c'e' ^^^^ '^^ ^^"^ ^0 being the ground-state annihilation and creation op- 
erators, as suggested in Ref. [H] based on results in [19], and which has been utilized 
by others, see, e.g., [20]. The "law of large numbers" considered by us in order to make 



the constant to vanish in the Bose-condensed phase has the form (6N /Ny j oc N 
see, e.g., ^ Eqs. (3.55) & (3.57)]. 



We substitute x'^li ^) = ""-^ J2k Ck [S {Tiu — Ak(q)) — 5 {huj + Ak(q))] from [T7] for the 
IQG with Ak(q) = £^|k+q| — £^fc, and obtain for q 7^ , 



(g>0) 



(3) 



where Ck denotes the thermal-average fraction of particles having momentum ^k, 



9s 



9s 



(4) 



Here 7] = +1,-1,0 refer to BE, FD, MB gases, respectively. gs=2s + 1 is the spin- 
degeneracy factor for spin s, A = e'^^ is the fugacity, and the function (u{x) denotes the 
polylogarithm [121 US] of order u. The solution of the equation ji = jir, {n, T) =0 gives the 

characteristic temperature which can be expressed as T( 



with = h'^k"^/ (2m) and fcu = 2 {Gn^n/ gs)^^^ serving as units of energy and wave number, 
respectively. 

For further analytical discussions, Eq. (jS]) will now be recast into an appropriate form by 
(i) substituting k — )■ k + q and to get 2S{q), (ii) expressing exponentials in coth-functions 
in accordance with the first equation in (j4]), and (iii) using J2k {pk + C'|k+q|) = 2. The 
procedure finally yields a new form: 



'u/kB (6v/7?C3/2(^)/^ 



-2/3 



S{q) - 1 = !^ ^ C, qk+q| , (g > 0) 
9s k 



(5) 



from which we read 



1 



^Ee^''-[%)-l] = ^^2CoEC,e^''- + 



(6) 



with Co=Nq{T)/N denoting the fraction of particles which occupy the zero-momentum 

state. Introducing the thermal de Broglie wavelength A=y^27r/i^/3/m and the dimension- 
less function 



Fir) 



E 

k 



e^^-^Cfc = Co(T) 
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(7) 



which is related to the one-body density matrix [21j by n^^\r, r') = {n/gs) F{\r — r'|) , the 
condensed fraction Co{T) can be extracted from the normalization condition -F(O) = 1. 
The evaluation of the ^-integration thus leads to Co(T) = 5^ i9(Tc — T) 1 — (T/Tc)^^^ , 
with G(x) the Heaviside unit step and 6ij the Kronecker delta, in conformity with [17] 
and the condensed-IB G result [H] Chap. 3.2]. 

Equations ([5]) and constitute to be our basic results which are valid at all temper- 
atures and for all ideal gases. Inserting Eq. ([6]) into Eq. ([T]), we find 



9[r) 



1 + 



F\r) - Cl{T) 



(8) 



which, in conjunction with Eq. ([7]), yields an expression in agreement with that discussed 
in Problem 4 of [51 §117] for an IBG at T < T^. Thus Eq. (jS]), which is valid at all 



temperatures for all IQGs, generalizes |6i Eq. (117.8)] whose validity is for an FD gas at 
all T but for a BE gas at T>Tc only. 

It seems pertinent to mention that for bosons g{r) is not simply the sum of condensate 
and non-condensate (or thermal) contributions as it is for the one-body density matrix, 
Eq.(l7]), or the density-response function jXTjEl]- The presence of the Fourier convolution 
in Eq.(l5]) has resulted into Eq. ([8]) wherein the thermal contribution [F{r) — Co(T)] 
appears as a factor in the second term on the right hand side. This factor vanishes in the 
limit T — )■ and hence is responsible for the "flattening", g{r) 1, observed at T <^ T^. 

There is another significant aspect regarding the derivation of our results in context of 
the GCE used here. On the dynamic route leading to Eq. ([3]), and thereby Eq. ([5]), we 
neither needed to nor did we use the Bogoliubov prescription which replaces and 
by c-numbers. For a consistency check, we took recourse to the static route (not elabo- 
rated here) starting from S{q)= SN^SN^ ^ /N with iVq= X^k.o- ^.ak+q,o-. We obtain for 
a uniform fluid the counterpart of Eq. (j5]) wherein the right hand side contains additional 
terms representing correlations of number fluctuations, 6 (o^o-'^k'.cr') ^ ('^k'+q o-''^k+q,cr) )• 
However, for any q 7^ 0, these extra terms vanish for ideal gases, irrespective of population 
of any single-particle state. Hence it is comforting to note that static and dynamic routes 
lead to exactly the same result. 

Analytic Expressions for F{r): For a BE or an FD gas, the integral in Eq. ([7]) can be 
carried out analytically in the region — 00 < yU < 0, i.e. for < A < 1. This region covers the 
complete domain of the IBG while it describes only the high-T domain (Tq <T < 00) 
of the ideal Fermi gas (IFG). On series expansion of Ci(^) ^ind subsequent term-by-term 
integration, we get 



An alternative form, equivalent to Eq. ([9]) and most suitable for evaluation at r^A, is 



which results from series expansion of the exponential function in Eq. ([9]) and subsequent 
interchange of summations. Also, we have used the relation obtained on implementing 
F(0)=1 in Eq. (E]). 

For analytical evaluation of the integral in Eq.([7]) for 1]= — 1 in the region where /i>0, 
i.e. for A > 1, which describes the low-T domain (0<T<To^^'') of an IFG, we split the 
integral into a sum of two integrals over intervals (0, and (^/5/I, 00), respectively. 

In the latter integral, the polylogarithm Ci(-2) with z = —e'^^"'^ can be expanded into a 
power series in z, since < < 1. In the former integral, where \z\ > 1, we apply the 
identity (i{z) = (i{l/z) — \n{—z) vahd for z^ (0, 1) and subsequently expand ^1(1/2;) into 
a power series in 1/z, since \l/z\ = l/\z\ < 1. Term-by-term integration of the resulting 
infinite sum finally yields for low-T IFG (/U>0) : 
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where ji{x) is the spherical Bessel function of the first kind and order 1, and erfc(x) is 



the complementary error function. Also, kp = h 2m/i_i(n, T) denotes a generalized 

Fermi wave number with limT^o kp = kF = {Gir'^n/ gsY^^ and is a measure of the chemical 
potential. In the low-temperature limit, i.e. for T— J-O and A— )-cxd, one easily retrieves 
the following from Eq. f lTTj) : F™(r) 3 ji(fcF'")/(/i^F'"), the expression given, see e.g., 

in Ref. flOl. 



Pair Distribution Functions: From Eq. (|8]), one readily finds for a dilute quantum gas, 
g^^{r) = 1 which coincides with the classical ideal-gas result. One also deduces the results 
5^(0) = ! + {ri/gs)[l — Cq{T)] and g{oo) = l leading to the following bounds: 
^^^(0) = 1 - l/gs<g^^{r)<g^^{oo) = l and ^^^(oo) = l<^^^(r)<^^^(0) with 



l + l/gs if T>T, 

if T <Tr 



l + {l/gs) \2 {T/T,f'-{T/Tf 



:i2) 



The small- and large-r asymptotic behaviors for T <Tc are obtained as g^^{r — 0) = 
g^^iO) - (2vrrV<7,A2) (T/T,)^/^ [C(5/2)/C(3/2)] + O ((r/A)^) and g^^{r oo) = 1 + 
2Co(T)/ (nA^ r)+gs/ (n^A^ r^). We find that the latter asymptote improves the expression 
given in [9l Eq.(21)], and is in agreement with [6l p. 359]; the expression given in [9] would 
be valid only at T^T^ whereas the validity of ours is in the entire range 0<T<Tc. Also, 
for T^O, one gets = - 1)/^, + (fc^rVS^,) [1 - 3A;2 r735 + ^(r^)] wherein 

the first two terms on the right-hand side give the result as obtained by Lee and Long 
[To] while discussing the static structure for an ideal electron gas at T=0. And for large 
distance r, we get g^^{r — )■ oo) = 1 — 9 [cos'^{kFr) — sm{2kFr) / (kpr)] /{ggkpr^) + 0{r~^) 
approaching unity as r~^ by damped oscillations, which improves on a result given in [6l 
§117]. 

On substituting in Eq. the high-T asymptote of F(r) deduced from Eq.dH]), the 
Gaussian form, 

g{r) 1 + — expl — I (13) 

is obtained. In fact, the asymptote ( 1T3|1 generalizes to the spinor gases the earlier results 
derived by others for zero-spin particles, see e.g., [22]. Although the experimental results 
for a thermal bosonic gas were fitted by Schellekens et al. [5] using an expression like Eq. 
(fT3|) . a detailed analysis together with the theoretical plots showing contrasting behavior 
for spinor bosons and fermions is lacking. 

The computed values of (2s + 1) [g{r) — 1] = rj [F^(r) — CliT)] are depicted in Fig. 
□ for a set of six reduced temperatures T* = fcBT/e„=0.491, 0.164, 0.114,0.104, 0.055 
and 0.011 corresponding to T/Tc=4.5, 1.5, 1.05, 0.95, 0.5, 0.1, respectively. The value 
g^'^^{r)=l of an ideal dilute quantum gas coincides with the abscissa at all T and, it 
is to be noted, that the pair-correlation properties of an IQG differ qualitatively from 
those of the corresponding dilute gas, even at high temperatures. As displayed in FigUl 
g^^{r) > 1 and g™{r) < 1 at all T. We term the appearance of a bump (dip) as the 
"Bose pile" ( "Fermi hole" ) which reflects the statistical "effective" attraction (repulsion) 
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Figure 1: (Color online). Pair correlation functions of ideal BE (above ab- 
scissa), FD (below abscissa), and MB (coinciding with the abscissa for all 
r) gases at reduced temperatures T*= 0.491 (a, red), 0.164 (b, green), 0.114 
(c, blue), 0.104 (d, short-dashed), 0.055 (e, dot-dashed), and 0.011 (f, long- 
dashed). T=0 results: bosons (coinciding with abscissa), fermions (lowermost 
curve, black, shadowing the long-dashed line). 



in ideal gases, an interaction which weakens with increasing spin. For 1, one would 
recover the result g{r)=l true for a classical ideal gas. 

It is found in low-r region that g^^{r) monotonically decreases (increases in magnitude) 
with increase in r as T decreases. However, for bosons in that domain, g^^{r) increases 
as T decreases for T > Tc and the trend reverses with decrease in T below T^. Let us look 
further into the behavior of fermionic curves at all T and bosonic at T > Tc, i.e. those 
marked (a) to (c) above abscissa. It can be seen that the pile (hole) being generated by 
rotating the Gaussian curves about the ordinate merely narrows down retaining its original 
spin-dependent height (depth) as T is increased. They ultimately acquire congruent 
bell-shaped forms (cf. Eq.(TT3l)) in the high-T regime having width (at half maximum) 



= Yln2/27r A ^ 0.332 A. Also, it can be seen that g^^ir) and g^^{r) show different 
behaviors at T = 4.5 Tc as opposed to studies wherein it was found that response functions 
|17j . dynamical structure factors [23], and momentum distribution functions [T7] of BE, 
FD, and MB gases are essentially independent of statistics at this temperature. However, 
in the limit T— j-oo, we get g{r)=l independent oi statistics, as expected. 

The bosonic curve marked (c) clearly demonstrates that g^^{r) becomes long-ranged 
as T reaches in the close vicinity of Tc from high-T side. However, unlike this, g^^{r) is 
of much shorter range at all T. It is estimated that the correlation length C,fd{T) < 10 
which presents a measure of the largest distance at which fermion pairs are still correlated. 
If T is raised, ^fd(7') further decreases and at high T it is of the order of A (cf. Eq. (fT3|) ). 

It is tempting to compare both the curves marked (a) in Fig. [1] with Fig. 2 of Ref. 
j3] wherein the HBT effect showing bunching and antibunching for ultracold atoms of 
^He* and ^He* have been depicted. The striking resemblance between the theoretical 
and experimental plots is remarkable. The direct comparison of our curves with the 
experimental ones will have to await availability of data free from the uncertainty and 
systematic errors in measurements, mentioned by Jeltes et al. [5]- However, the interesting 




physics is unfolded in the BEC phase represented by the curves (d) to (f). The height of 
the Bose pile goes on decreasing and the curve gets increasingly flattened as T — (cf. 
Eq. (fT2|) ). The curve (f) for T = 0.1Tc is almost flat and ultimately the plot for an IBG 
at T = coincides with the abscissa analogous to the Bose-condensed phase experimental 
result IM] revealing that the system is completely coherent. In fact, it is just like the 
situation in a single-mode laser in which the photons are not bunched [25]. Furthermore, 
our studies depict temperature-dependent aspects of bunching and antibunching. 
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